Tensor Microwave Anisotropics from a Stochastic Magnetic Field 
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We derive an expression for the angular power spectrum of cosmic microwave background 
anisotropies due to gravity waves generated by a stochastic magnetic field and compare the result 
with current observations; we take into account the non-linear nature of the stress energy tensor of 
the magnetic field. For almost scale invariant spectra, the amplitude of the magnetic field at galactic 
scales is constrained to be of order 10 -9 Gauss. If we assume that the magnetic field is damped 
below the Alfven damping scale, we find that its amplitude at 0.1h _1 Mpc, B\, is constrained to be 
B x < 7.9 x 10" 6 e 3n Gauss, for n < -3/2, and B\ < 9.5 x 10" 8 e°- 37n Gauss, for n > -3/2, where 
n is the spectral index of the magnetic field and Ho = 100/ikm s _1 Mpc _1 is the Hubble constant 
today. 
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PACS Numbers : 98.80. Cq, 98.70.Vc, 98.80.Hw 
I. INTRODUCTION 

The past few years have seen a tremendous surge of 
interest in the origin and evolution of galactic magnetic 
fields jj] . A number of mechanisms have been proposed 
for the origin of the seed fields, ranging from inflationary 
mechanisms M, cosmological phase transitions [|| to as- 
trophysical processes Q]. Much progress has been made 
in trying to disentangle the various non-linear processes 
which may be responsible for the growth of such a seed 
field in the very early universe in particular the interplay 
between the magnetic field and the primordial plasma 
HD an d the importance of turbulence . 

Given a small seed field at late times, two different 
mechanisms can cause its amplification to magnetic fields 
of order 10 -6 Gauss observed in galaxies: adiabatic com- 
pression of magnetic flux lines can amplify a seed field 
of order 10 -9 Gauss to the present, observable values; 
the far more efficient (and controversial) galactic dynamo 
mechanism may be able to amplify seed fields as small 
as lCr 20 Gauss [| or even 10 30 Gauss in universe with 
low mass density ||. Clearly, to make some progress in 
identifying which one of these mechanisms is responsible 
for galactic magnetic fields, one would like to find a con- 
straint for the seed field before it has been processed by 
local, galactic dynamics. 

The obvious observable for such a constraint is the 
cosmic microwave background (CMB). It is interesting 
to note that a field strength of 10 _8 Gauss provides an 
energy density of Qb = B 2 /{%-Kp c ) ~ 10 _5 f2 7 , where 
f2 7 is the density parameter in photons. We naively ex- 
pect a magnetic field of this amplitude to induce per- 
turbations in the CMB on the order of 10~ 5 , which are 
just on the level of the observed CMB anisotropies. It is 
thus justified to wonder to what extent the isotropy of 
the CMB may constrain primordial magnetic fields. Our 
order of magnitude estimate makes clear that we shall 
never be able to constrain tiny seed fields on the order of 



10 _13 Gauss or less in this way, but primordial fields of 
10 _9 Gauss may have left their traces in the CMB. 

A number of methods have been proposed in the past 
few years for measuring a cosmological magnetic field us- 
ing the CMB: the effect on the acoustic peaks || , Faraday 
rotation on small |lo) and large [|ll| scales and vorticity 
Jl2| , ^3t can all lead to observable anisotropies in the CMB 
if the primordial magnetic field strength is of the order of 
1CP 9 to 1CP 8 Gauss. The most stringent bound from the 
CMB presented thus far was for the case of a homoge- 
neous magnetic field |l4j] ; the authors use the COBE data 
to find the constraint B < 6.8 x 10 _9 (il /i 2 ) 1 ^ 2 Gauss 
where the Hubble constant is Hq — 100/ikm s~ 1 Mpc _1 
and Qq is the energy density in units of the critical value. 
Although there is no fundamental reason to discard the 
possibility of a homogeneous magnetic field, all physi- 
cal mechanisms proposed to date lead to the presence 
of stochastic magnetic fields with no homogeneous term; 
in this paper we consider such fields. For these types 
of configuration one is allowed to have fluctuations on 
a wide range of scales and the magnetic field will serve 
as a non-linear driving force to the metric fluctuations; 
in the parlance of cosmological perturbation theory, the 
magnetic field evolves as a stiff source, without being 
affected by the fluid perturbations (back reaction) 
which may be induced. 

Stochastic magnetic fields have also been considered 
in p2| , where the CMB anisotropies due to the induced 
fluid vorticity has been analyzed. Here we determine 
gravitational effects of the magnetic field. For simplicity, 
and to allow for a purely analytical analysis, we constrain 
ourselves to tensor perturbations. Similar contributions 
are also expected from vector and scalar perturbations 
which then would add to the final result. In this sense 
the anisotropies computed here are a strict lower bound 
(underestimating the true effect probably by about a fac- 
tor of three). 

The main result of this work is that one can obtain 
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reasonably tight constrains for scale invariant magnetic 
fields; For causally generated magnetic fields the con- 
straints are weaker and are strongly dependent on the 
evolution of the magnetic field in the radiation era on 
small scales. 

For simplicity we concentrate on the case flo = 1. 
Througout, we use conformal time which we denote by 
77. Greek indices run from to 3, Latin ones from 1 to 3. 
We denote spatial (3d) vectors with bold face symbols. 
The value of the scale factor today is 0.(770) = 1. 

II. THE STRESS TENSOR OF THE MAGNETIC 
FIELD 

During the evolution of the universe, the conductiv- 
ity of the inter galactic medium is effectively infinite. In 
this regime we can decouple the time evolution from the 
spatial structure: B scales like B(t7, x) = Bo(x)/a 2 on 
sufficiently large scales. On smaller scales the interac- 
tion of the magnetic field with the cosmic plasma be- 
comes important leading mainly to two effects: on in- 
termediate scales, the oscillates like cos(w Akrj), where 
va = B 2 /(4:TT(p + p)) 1 / 2 is the Alfven velocity and on 
small scales, the field is exponentially damped due to 
shear viscosity || . 

We will model Bq(x) as a statistically homogeneous 
and isotropic random field. The transversal nature of B 
leads us to 

(B z (k)B*(q)) = S 3 (k - q)(<^ - k l k J )B 2 (k) . (1) 
where we use the Fourier transform conventions 
Bj(k) = J d 3 xexp(ix ■ k)_Boj (x) , 
B o,(x) = ^3 J d 3 fcexp(-zx.k)S,(k) . 

The Alfven oscillations modulate the initial power spec- 
trum by a factor 

B 2 (k) -> B 2 (k) cos 2 (v A kr]) . 

This can be approximated by a reduction of a factor 2 
in the power spectrum on scales with VAkr/ ~ 1. But 
as we shall see, our most stringent constraints will come 
either from very small scales where the spectrum is ex- 
ponentially damped or from much larger scales where 
oscillations can be ignored. We will incorporate the ex- 
ponential damping by a cutoff in the power spectrum at 
the damping scale. 

Let us investigate the consequence of causality for the 
spectrum B 2 (k). If B is generated by some causal mech- 
anism, it is uncorrelated on super horizon scales, 

(B i (x,T7)B i (x / ,» 7 )> = for |x-x'|>2?7. (2) 



Here it is important that the universe is in a stage of stan- 
dard Friedmann expansion, so that the causal horizon 
size is about r\. During an inflationay phase the causal 
horizon diverges and our subsequent argument does not 
apply. In this somewhat misleading sense, one calles in- 
flationary perturbations 'a-causal'. 

According to Eq. (||), (Bi(x, r])Bj(x', rf)) is a function 
with compact support and hence its Fourier transform is 
analytic. The function 

(Bi(k)B;(k)) = (6ij - kk 3 )B 2 (k) (3) 

is analytic in k. If we in addition assume that B 2 (k) 
can be approximated by a simple power law, we must 
conclude that B 2 (k) oc k n , where n > 2 is a even integer. 
(A white noise spectrum, n = does not work because 
of the transversality condition which has led to the non- 
analytic pre- factor Sij — kikj.) By causality, there can be 
no deviations of this law on scales larger than the horizon 
size at formation, r?j n . 

We assume that the probability distribution function 
of Bo is Gaussian; although this is not the most general 
random field, it greatly simplifies calculations and gives 
us a good idea of what to expect in a more general case. 

The anisotropic stresses induced are given by the con- 
volution of the magnetic field, 

T iP( k ) = ^ / AA(q)S,*(k-q) 

-^(q)Bf(k-q)«% . (4) 

With the use of the projection operator, Pij — Sij — kikj 
we can extract the tensor component of Eq. (Q) , 

11? ] = (P?P" ~ (l/2)P l3 P ab )r ab , (5) 

tracelessness, orthogonality and symmetry force the cor- 
relation function to be of the form 

<n<f (M)n^*(k',*)) = |n B (fc,<)| 2 M y7m <5(k-k') 
(r4f (M)njf *(k',*)> = 4|iL3(M)l 2 <5(k-k'), (6) 

were we make use of the tensor basis, M.: The corre- 
lator on an isotropic tensor component has always the 
following tensorial structure, 

M-ijlm = 5u8j m + Si m 6jl — SijSl m + fc 2 (Sijkik m + 

^imkikj 3ukjk m Si m k[kj djikikjji 
fijmkiki) + k kikjk\k m . (7) 

We now determine the function |IlB(A:,i)| 2 in terms of 
the magnetic field. Using Wick's theorem we have 

(s l (k)s;( q )s, i (s)i3;„(p)) = 

(B l (k)B*( q ))(B n (s)B; n (p)) + 
(B i (k)B* n (s))(B j ( q )B* m (p)) + 
(Bi(k)B* m (p))(B n (s)B*( q )) . (8) 
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The problem reduces itself to calculating self convolu- 
tions of the magnetic field. The power spectrum of 
Eq. (@) is 

(7"«(k,r/)r^*(k',?7)) = 

J cfqj d 3 p(B l (q)i? J ((k-q)) J B i (p)i? m ((p-k'))> - 

<5(k - k') J d 3 q B 2 (q)B 2 (\k - q|) X 

[(Su - q t qi)(S jm - (k - q) J (k - q) m ) + 

(S im -MmWjl ~ (k- q)j(k- q)/)] ■ (9) 

Using Eqs. (§J|) and ([§, this leads to |n B | 2 = f(k) 2 /a s , 
where 

/(fc)2 = (8^F / A^foWlk - + 27 2 + 7 2 /3 2 ) , 

(10) 

with 7 = k • q and (3 — k ■ k — q. 

It remains to define B(k) from Eq. (Q). We shall pa- 
rameterize it in terms of an amplitude and a scale depen- 
dence through 



B 2 (k) 



'ygj B^" forfc<fc c 







otherwise 
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n > —3/2, the spectrum of the energy momentum ten- 
sor becomes white noise, independent of n. Only the 
amplitude which is proportional to (Xk c ) 2n depends on 
the s pec tral index. This is due to the fact that the inte- 
gral (10) is dominated by the contributions at very small 
scales, k c 3> k. The induced Cg spectrum from gravity 
wave will therefore be independent of n for n > —3/2, 
and obey the well known behavior Ci oc t of a white noise 
source. 

To simplify, we just consider the dominant term and, 
in order not to artificially produce a singularity at n — 
—3/2, we drop the factor l/(2n -I- 3). Given the intent of 
this paper (to constrain the amplitude of the magnetic 
field) we will include a factor of 10 _1 in our final re- 
sult, guaranteeing that we are not overestimating CMB 
anisotropics. The singularity at n = —3 is real. It is 
the usual logarithmic singularity of the scale invariant 
spectrum. 



III. THE CMB ANISOTROPIES 

Armed with the structure and evolution of the stochas- 
tic magnetic field we can now proceed to calculate its ef- 
fects on tensor CMB anisotropies. The metric element of 
the perturbed Friedman universe is given by 

ds 2 = a 2 (-q)(-dri 2 + (8. tj + 2h ij )dx i dx j ) , 



The normalization is such that (Bq(x)Bq(x))\x = B 2 
where the quantity in brackets represents the averaged 
magnetic field smoothed over a comoving length scale A. 
Note that we have assumed that the cutoff scale today is 
smaller than A. 

We require n > —3 so as not to over-produce long 
range coherent fields; we shall see that for n = — 3 we 
obtain a scale invariant spectrum of CMB anisotropies. 

We have included a short wavelength cutoff to take into 
account the exponential damping due to shear viscosity 
in the cosmic plasma || . The mean energy density due to 
such a magnetic field, which is an appropriately weighted 
integral of Eq. (|TT|) , will be strongly dependent on the cut 
off when n > —3. 

Using Eqs. (|ll| ) and (10) we can calculate /. The inte- 
gral cannot be computed analytically, but the following 
result is a good approximation for all wave numbers k 



f(k) 



(2nf 



A 2„+6 B 4 



16 r 2 m±}(2n + 3) 



2n+3 



(12) 



This result seems to have a singularity at n = —3/2 
which is however removable. The first term dominates 
if n > —3/2 and while second term dominates if the 
opposite inequality is satisfied. For n > —3/2, the grav- 
ity wave source is therefore white noise and its ampli- 
tude is determined by the upper cutoff, k c . Note that if 



where h\ — and h\k % = for tensor perturbations. The 
magnetic field will source the evolution equation for 
through 



hij + 2-hij + k 2 hij = 8TrGU\f ] 



(13) 



Such a gravity wave induces temperature fluctuations in 
the CMB due to the fact that the photons move along 
the perturbed geodesies |15 



AT f Vo ■ - ■ 

-(?7o,x,n)=/ /i l:) (x(?7),?7)nVcZf7 . (14) 



T 



Here 77* denotes the (conformal) time of decoupling of 
matter and radiation due to recombination. We want 
to compute the angular power spectrum of ^r, the Ci, 
defined by 



AT AT 

— (") — (n'))n-n'= M r 



P^(2^ + 1)C,P,(aO 



The CiS are solely determined by the power spectrum of 
metric fluctuations. Defining 



(h^'ik'^h^iKr,)) = \H(k, v )\ 2 M ljlm 6(k-k>) 
one can derive a closed form expression for Cg (see M): 
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c t = 4^4 / dkk 2 \I(£, k)\ 2 £(£ - 1){£ + 1)(£ + 2) , (15) 



{Km - vW 



(16) 



where jt denotes the spherical Bessel function of order £. 
We solve equation (|13|) using the Wronskian method; in 
terms of the dimensionless variable x — kr\. The homoge- 
neous solutions are the spherical Bessel functions jo , yo 
in the radiation dominated era, and j± jx , y\ jx in the 
matter dominated era respectively. We assume that the 
magnetic fields were created in the radiation dominated 
epoch, at redshift z\ n . We then match the general in- 
homogeneous solutions of Eq. (|l^) at the time of equal 
matter and radiation, r\ eq . Due to the rapid falloff of the 
source term in the matter dominated era, the perturba- 
tions created after r] eq are sub-dominant, and we find for 
the dominant contribution at r\ > 77, 



eq 



H{k, t) ~ 4irGn 2 z eq In ( — ) kf{k) J -^- . (17) 



c e q j kfj 

Inserting this result in Eq. (EBh, we obtain 



I = AirGr/QZeq \n(z in / z eq )f(k) / dx 



J 2(x) jljxp - X) 
X (xq — x) 2 ' 

(18) 



where x = krj, x* = kr/* and xq = kr]Q. For wave numbers 
which are super-horizon at decoupling, x» < 7r, the lower 
boundary in Eq. ( |l8| ) can be set to 0. The remaining 
integral cannot be expressed in closed form, but is well 
approximated by |p0[: 



Xo , ji{x Q - x) 
dx 7 ^2" 

X (Xq — X) z 



7T f X " Jb/2{x) J e+ 1/ 2 (XQ ~ X) 

2 J X x 3 / 2 (x ~x) 5 / 2 ~ 
0.7% r X0 dx J 5/2(x) %i/ 2 (io - x) 
2 Jo x (x - x) 3 



7ttV~£ t , „ 
— -g-Je+aixo)) . 



(19) 



The third integral above can be expressed in closed form 
( |pl| , number 6.581.2), and is reasonably well approx- 
imated by the last expression, we have checked the ap- 
proximation numerically for I < 200 and varying xo- 

We can now do the integrals in Eq. ( [l5| ) analytically 
to obtain 



fQ ~ A 



\\ 2n + 6 2 



170 



Sit 



(k c r) 



,271+3 /j3 



for n > —3/2, and 



2 Ci ~ A 



2n+6 



T[l - 2n] 



(n + 3) r 2 [l -n]2(!- 2 ") 



6+2n 



(20) 
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for —3 < n < —3/2, where 
A = 5 x lti- i {2irfz 2 e In 2 



BtGW 

r 2 



■ n+3 1 



3 x 10" 



10- 9 Gauss 



h\ 2 (z in /z eq ) 



r 2[«±3] 



(22) 



IV. RESULTS 

Eqs. (^) and (|2^) are our main result. They allow us 
to limit a possible primordial magnetic field by requiring 
it not to over produce fluctuations in the CMB. Since 
the fluctuations induced grow with £ for all values of the 
spectral index —3 < n, we obtain the best limits for 
large values of t. We shall be conservative and assume 
an upper bound of l 2 Cg |£=5o< 8.5 x 10~ 9 Q. Given 
that we are interested in galactic and cluster scales we 
fix A = 0.1 h _1 Mpc for the remainder of this paper. In 
Fig. 1 we show the limit on a stochastic magnetic field 
as a function of the spectral index n, using the damping 
scale given below as cutoff. 

We now focus on a few particular cases of interest and 
in doing so we will derive an analytic expression which ap- 
proximates the upper bound of B\ over the whole range 
of n. 

Scale invariant magnetic field: From Eq (|2l]) we see 
that the result is independent of the cutoff. In the limit 
where n — ► — 3 we find that 



10 J Gauss 



(23) 



1. e. of the same level as other constraints ]9j-|li 

Causal magnetic field: For this scenario we have, as 
explained above, n > 2; we shall consider the case of n = 

2. For instructive purposes let us first consider a k c which 
is independent of the magnetic field. The constraint is 
then 



B x £ In" 



(kcVo) 4 Gauss 



(24) 



The cutoff k c will depend on the plasma properties and 
evolution; even though the conductivity a of the cosmic 
plasma is very large, it is nevertheless finite. One actually 
finds |l9| that a = aT, where the parameter 1 < a < 7 
is slowly temperature dependent. By Ohm's law, mag- 
netic fields on small enough scales are exponentially sup- 
pressed, B oc exp (— ak 2 r]/4:Tra) , leading to a damping 

scale, kd(v) = (^cra/r]) 1 / 2 = (77 x 2 x lO^cm) -1 / 2 . This 
scale is smaller than the comoving horizon scale for all 
temperatures below the Planck scale. On scales smaller 
than \/kd{rj eq ), the induced gravity waves have damped 
away even before matter and radiation equality. Since the 
sourcing of gravity waves after equality is negligible, the 
damping scale relevant in our problem is k c = kd{rj eq ), 
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2 x 10 13 /i 2 Mpc 



(25) 



If we insert this this damping scale in Eq. ( p24| ) , we obtain 

B x ~ 10- 29 Gauss. 

A more realistic scenario is to assume that the mag- 
netic field will be damped by electron viscosity. To 
proceed with the analysis we shall split the stochas- 
tic magnetic field into a high-frequency component and 
a low- frequency component; the scale which separates 
the two is the Alfven scale at equality A ,4 = VA^eq 
where Va is the Alfven velocity, Uf = {B 2 )/{Att{p + p)). 
From equation 4 of |l^] we see that the inhomogeneous 
magnetic field will obey a damped harmonic oscillator 
equation, with a time dependent damping coefficient, 
D = 0.2fc 2 Z 7 (l + z) (l^ is the physical photon mean free 
path) and frequency cj(r?) = V\k 2 - (D/2) - {D/2f. 
Within this setting we can estimate the damping scale of 
the magnetic field in the oscillatory regime of this sys- 
tem; the amplitude of the effective homogeneous mag- 
netic field, Ba , which is responsible for the Alfven waves 
is related to B\ through 



100 = 



3 

d 
o 



m 




FIG. 1. The upper bound B\ as a function of spectral in- 
dex, n. We assume z in / z eq — 10 9 and A = 0.1h~ 1 Mpc 



B\^B\ 



n+3 



which leads to 



B A 



Bx 



10- 9 Gauss 



A10~ 9 Gauss 
3.8 x IO-^BaVc 



(13ft -1 )Sfu(r 9 Gaiiss (26) 



We shall define the damping scale to be the scale at which 
one e-fold of damping has occured by equality. From 
Io' q %dr] = \ one finds 



k c = 4.5Mpc~ 



(27) 



For this estimate to be valid, the system must be in the 
damped oscillatory regime (as opposed to overdamping 
regime), i.e. u>{rj eq ) > 0; this condition is satisfied if 
Ba > 5.5h~ 2 x 10~ 9 Gauss. We find that indeed this is 
the case in the range of interest. 

Combining Equations (|l|) and @ and as- 

suming a formation redshift of = 10 15 (although the 
final result is very weakly dependent z.;„) we find that an 
an approximation to the bound is 



B x < 7.9 x l(r 6 e 2 ' 99n Gauss, 
B x < 9.5 x l(r 8 e a37n Gauss, 



for n < -3/2, 

for n> -3/2, (28) 



The upper bounds corresponding to Eq. ^8| represent a 
reasonable fit to Figure 1. As one can see, the constraint 
on a causal magnetic field is well above 1CP 9 Gauss. 

Throughout this derivation we have assumed that we 
can estimate the damping scale of the magnetic field by 
looking solely at the Alfven modes. A linear analysis of 
the remaining degrees of freedom also indicate that the 
magnetic field will be damped at the same scale as in 
Equation (|27j). It is possible that non- linear effects may 



prevent the tangled magnetic field from damping at this 
scale but an accurate quantitative analysis is still lacking. 

Inflationary magnetic fields: Broken conformal invari- 
ance allied with the an inflationary period will create 
large scale magnetic fields. Bertolami and Mota [|| es- 
timate the spectral index in such a mechanism to lie in 
the range around 0. We are then clearly in the regime 
where the cutoff is important. Using the Alfven damping 
scale at equality and assuming magnetic field generation 

at 10 10 GeV, we find B x ~ (1CT 7 - l(T 8 )Gauss for n vary- 
ing from —0.5 to 0.5. A similar result can be obtained 
for the model of Gasperini et al M. 



V. DISCUSSION 

Our calculation differs from most of the recent work on 
the impact of primordial magnetic fields on structure for- 
mation: In estimating the CMB anisotropies we do not 
split the magnetic field into a 'large' homogenous mode 
and a 'small' fluctuation. The magnetic field then affects 
metric perturbations quadratically. This has two effects. 
Firstly it allows us to consider the magnetic field as a 
stiff source, and discard (within the MHD approxima- 
tion) the backreaction of the perturbations in the cosmo- 
logical fluid. Indeed if we were to consider backreaction 
then we would know a priori that we would be generat- 
ing unacceptable perturbations in the cosmological fluid. 
Another way of phrasing this is that the magnetic field 
itself is 1/2 order perturbation theory, while its energy 
momentum tensor and consequently the induced metric 
perturbations are first order perturbations. The MHD 
backreactions on B would be 3/2 order and may thus be 
neglected in linear perturbation theory. We point out, 
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however that, to obtain an estimate of the damping scale 
due to the viscosity in the MHD we had to consider a split 
between long wave length and short wavelength fluctua- 
tions in B. 

Secondly, the stress energy tensor being quadratic in 
the magnetic field, leads to a 'sweeping' of modes: large 
wavelength modes in T M „ will in general be affected by 
all scales of the spectrum of B As we have seen 

in the causal case, the small wavelength behaviour of 
the magnetic field totally dominates the large wavelength 
pertubations. In || the magnetic field is modeled as 
B = B + B^'(x) where B is a homogeneous term; the 
stress energy tensor is then given by terms of the form 
BiBf\ which are linear in the stochastic component. 

A few comments are in order with regards to our result. 
Note that we are considering a specific class of models, 
where the magnetic field seed is created at some well 
defined moment in the early universe and then evolves 
according to the MHD equations. If the magnetic field 
is being constantly sourced throughout the radiation era, 
then our calculation is not valid. An example of such a 
scenario was proposed by Vachaspati || where magentic 
fields are sourced by vortical imprints from an evolving 
network of cosmic strings; although the scaling behaviour 
of source may lead to B cx a~ 2 , the effective damping 
scale will be of order the horizon much larger than the 
Alfvcn damping scale. Another possibility has been put 
forward in |7]], where the onset of turbulence induces an 
amplification of power on large scales but a supression 
of power on small scales. This would further increase 
k c but the results are still too qualitative to be properly 
included in an analysis such as ours. 
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